Annual temperature, body size, and sexual size dimorphism in the evolution of Pyrgomorphidae

Abstract In many animal species, larger body size is positively correlated with male mating success and female fecundity. However, in the case of insects, in high seasonality environments, natural selection favors a faster maturation that decreases the risk of pre‐reproductive death. However, this advantageous adaptation comes at a tradeoff, resulting in a reduction in body size. Maturation time is influenced by environmental factors, such as temperature, season length, and food availability during the rains. The geographic variation in these parameters provides an opportunity to study their impact on the adaptive evolution of body size in Pyrgomorphidae grasshoppers. These grasshoppers exhibit remarkable variation in body size and wing development and can be found in diverse plant communities across Africa, Asia, Australia, and tropical America. In this study, we utilized a phylogenetic approach to examine the evolution of body size, considering climatic factors, and the influence of sexual selection on size differences between males and females. We found a positive correlation between mean annual temperature and sexual size dimorphism (SSD). Remarkably, species exhibiting a strong bias toward larger females were found to be adapted to regions with higher temperatures. In the Pyrgomorphidae family, an intermediate body size was identified as the ancestral trait. Additionally, winged male and female grasshoppers were observed to be larger than their wingless counterparts. Despite the potential conflicting pressures on body size in males and females, these grasshoppers adhere to Rench's Rule, suggesting that sexual selection on males' body size may explain the evolution of SSD.

In temperate regions or high altitudes, the survivorship of insects declines sharply at low temperatures, at the end of the fall and in the winter (Amarasekare & Sifuentes, 2012), and natural selection favors shorter development times to reduce the risk of pre-reproductive death, but this comes at the cost of smaller body size (Thornhill & Alcock, 1983;Dingle et al., 1990;Abrams et al., 1996;Ramírez-Delgado & Cueva del Castillo, 2022).Conversely, in regions with lower latitudes and altitudes, where temperature variations are minimal and plant growth is abundant, resulting in an abundance of food for herbivorous insects (Amarasekare & Sifuentes, 2012;Yom-Tov & Geffen, 2006), larger body sizes are favored by both natural and sexual selection in univoltine insect species (Kozlowski et al., 2004;Horne et al., 2015;Ramírez-Delgado & Cueva del Castillo, 2022).
In this study, using a phylogenetic framework, we analyzed the trends in the evolution of the body size of grasshoppers in the family Pyrgomorphidae.These grasshoppers are among the most recognizable orthopteran lineages in the world.Several species are known for their vibrant color patterns, elaborate sculptures, and large sizes (Figure 1) (Zahid et al., 2021).Nonetheless, most pyrgomorph species are actually cryptically colored.There is a large variation in their sizes, and there are winged and wingless species The species of the family are distributed in Africa, Asia, and tropical America (Mariño-Pérez & Song, 2019).Thus, they can be found in a wide variety of plant communities, including xerophytic, temperate, tropical deciduous, and rain forests.The variation in the niches they inhabit opens the opportunity to test the impact of the adaptive evolution on their body size.
We analyzed the evolutionary trends of body size among the pyrgomorphids, with special attention to the possible constraints imposed by the presence of wings.Furthermore, we explored the influence of climatic conditions on the evolutionary trajectory of both male and female body size.We also evaluated if the body size divergence in males and females can be attributed to sexual selection or fecundity selection testing the Rensch's rule.In several groups of animals, the magnitude of SSD increases with body size when males are larger than females, but it decreases with body size when females are larger than males (Rensch, 1950;Abouheif & Fairbairn, 1997;Fairbairn, 1997).This macro-evolutionary pattern is known as Rensch's rule (Rensch, 1950) and it has been associated with mating systems with high levels of mate competition (Abouheif & Fairbairn, 1997;Fairbairn, 1997).The converse trend, where female size varies more than male size, is less common, but seems to be the result of strong fecundity selection acting on females (Blanckenhorn et al., 2007;Foellmer & Moya-Larano, 2007;Webb & Freckleton, 2007).
Because environmental conditions can impact body size in grasshoppers, we expected that larger species of grasshoppers would be found in regions with low seasonality, high temperatures, and ample precipitation.These environmental factors contribute to the abundance of food for herbivorous insects, ultimately impacting their body size (Amarasekare & Sifuentes, 2012;Yom-Tov & Geffen, 2006).In addition, because resource allocation also plays a role, as flight can reduce resources for egg production, resulting in a tradeoff between flight and fecundity (Tigreros & Davidowitz, 2019), we hypothesized that winged males and females would exhibit smaller body sizes compared with their wingless counterparts.We also conducted allometric regressions of male body traits on female traits to assess whether sexual selection on male body size has been stronger than fecundity selection on females.As a result of stronger directional sexual selection on male body size, according to the allometric relationship predicted by Rensch's rule, a greater evolutionary divergence in male size than in female size would be expected (Berner & Blanckenhorn, 2006;Fairbairn, 2007;Stillwell et al., 2010).Thus, the slope of the allometric regression of male size on female size should be steeper than 1 (Rensch 1950;Fairbairn, 1997).

| Morphological traits
We measured museum specimens of taxa belonging to Pyrgomorphidae that we had previously collected and used to build morphological and molecular phylogenies deposited at the Texas A&M University Insect Collection (TAMUIC) (see Mariño-Pérez & Song, 2018, 2019;Zahid et al., 2021) del Castillo et al., 1999;Sanabria-Urbán et al., 2015).Considering the total number of individuals for each taxon, we averaged the values of the four morphological traits per species and sex (Table 1).The values obtained were then log-transformed prior to their utilization in the comparative analyses.Furthermore, we calculated the sexual size dimorphism index (SSDI) for each trait, following the methodology proposed by Lovich and Gibbons (1992) (Table 2).This index expresses SSD as [(length of larger sex/length of smaller sex) -1].For convention, the SSDi is arbitrarily changed to negative when males are the larger sex.Conversely, when females are the larger sex, the SSDi will be represented by a positive value.

| Climatic information
Because temperature and food availability can impact adult body size of grasshoppers, our study incorporated climatic parameters from the sampling localities of each specimen.Geographic coordinates were obtained using Google Earth 7.1.5(Google Inc 2014), while BioClim provided 14 climatic parameters related to temperature and precipitation (Hijmans et al., 2005).Due to the spatial resolution of the climatic data (30 s = ~1 km 2 ) and the scale of this study, we considered that the error associated to the sampling localities is negligible.Precipitation parameters were used as indicators of food availability, as supported by previous studies (Branson, 2010;Cueva del Castillo et al., 2015;Yom-Tov & Geffen, 2006).Our analysis focused on six climatic parameters directly associated with the growth and reproductive cycles of pyrgomorphids: the annual mean temperature, temperatures seasonality, the temperature of the wettest quarter of the year, annual precipitation, seasonal precipitation, and precipitation of wettest quarter of the year (Table 3).Due to the high correlation among climatic variables and the seasonal presence of nymphal and adult instar of these organisms, this approach allowed us to effectively capture the essential factors impacting their development and reprodcutive patterns.

| Phylogenetic reconstruction
To consider the phylogenetic effects in our comparative analysis (see below), we first reconstructed a phylogeny of Pyrgomorphidae based on mitochondrial genome sequences available on GenBank (Benson, 2017) that were generated by previous studies (see Table 4 for accession numbers and references of the used sequences).Our taxon sampling included three outgroups and 32 ingroup species.
We aligned all protein-coding genes individually considering the conservation of reading frames by first translating into amino acids in MUSCLE (Edgar, 2004) using default parameters in Geneious Prime 2023.1.Two ribosomal RNA genes (12S, 16S) were aligned in MAFFT using the E-INS-i setting, also in Geneious.Using SequenceMatrix (Vaidya et al., 2011) we concatenated all these individual alignments into a single matrix.We divided this matrix into a total of 15 data blocks (13 protein-coding and two ribosomal RNA).We performed a maximum likelihood (ML) analysis on the total evidence dataset (14,335 aligned bp and 44 taxa) applying the GTRCAT model to each partition and using RAxML 7.2.8 (Stamatakis et al., 2008) on XSEDE (Extreme Science and Engineering Discovery Environment, https://www.xsede.org)through CIPRES Science Gateway (Miller et al., 2011).Nodal support was evaluated using 1000 replications of rapid bootstrapping implemented in RAxML.
To estimate timing and rates of divergence in the phylogeny of Pyrgomorphidae, we performed a divergence time analysis using BEAST v.2.7.4 (Drummond et al., 2012).So far, there is no reliable fossil calibration point for the family.Therefore, we used the estimated age of the Pyrgomorphidae (109.35 mya) from Song et al. (2020) who used 11 fossil calibration points for inferring the divergence time estimate of Orthoptera.For our divergence time estimate analysis we used the best-fit models of nucleotide substitution and partitioning scheme recommended by PartitionFinder (Lanfear et al., 2016).Moreover, we relaxed the clock log-normal model for the clock model, the birth-death model with a uniform distribution as a tree prior, and a log-normal distribution as a distribution prior for the calibration point.We conducted two separate analyses each for 10 million generations, sampling every 5000 generations.To inspect convergence across independent runs, we used Tracer (Rambaut & Drummond, 2003-2009) and discarded 25% of each run as burn-in, and combined the two best trees that converged using LogCombiner (Rambaut & Drummond, 2002a).

| Body size ancestral state reconstruction
We conducted ancestral character state reconstructions to examine the evolutionary patterns of body size and SSD.We utilized the fastAnc function from the "Phytools" package in R 4.0.1 (R Core Team, 2015), as outlined by Revell (2013).This function estimates the ML value of a continuous trait for internal nodes and then F I G U R E 2 Anatomical characteristics measured in the grasshoppers studied species.
interpolates the states along the branches of the phylogenetic tree.
Additionally, we employed the contMap function in "Phytools" to visually represent the ancestral state reconstructions on the phylogeny for the traits under investigation.

| Climatic conditions and body size
We conducted PGLS regression models using R packages "nlme" (Pinheiro et al., 2019) and "caper" (Orme et al., 2013) to examine the influence of climatic variables on male and female body size and SSD.
We employed corBrownian and gls functions to convert phylogeny into correlation structure object and fit linear models, respectively.
TA B L E 1 Morphometric measurements of males and females of 32 species of the family Pyrgomorphidae.

| Body size differences between winged and wingless grasshoppers
We conducted phylogenetic ANOVA using the 'geiger' R package to explore potential differences between winged and wingless males and females.This package, as described by Pennell et al. (2014), incorporates a phylogenetic context by calculating the test statistic for ANOVA and simulating new sets of dependent variables on the phylogenetic tree.The null distribution of the test statistic is obtained under a Brownian-motion model, following the approach outlined by Garland et al. (1992).Males and females were analyzed independently.To estimate the phylogenetic signal of wing condition (winged or wingless) on the studied species we used the δ statistic (Borges et al., 2019), which is a real number that increases when a character evolves following the phylogeny and decreases when it evolves independently.For the estimation of the δ statistic we use the "ape" package (Paradis & Schliep, 2019) and the function delta with the default settings (Borges et al., 2019) in R. We test the statistical significance of δ statistic by generating a probability distribution of 200 random δ statistics (representing no phylogenetic signal) to compare the observed δ statistic.

| Rensch's rule
To assess whether male body size has diverged more than female body size in pyrgomorphid species (Rensch's rule), we employed the phylogenetic independent contrasts method (Felsenstein, 1985) via the "caper" R package (Orme et al., 2013) to account for species' phylogenetic non-independence (Harvey & Pagel, 1991).To ensure robust parameter estimates, we excluded outliers with studentized residuals exceeding ±3 (Jones & Purvis, 1997).Before proceeding with further analyses, we examined the independence of the standardized contrasts from their estimated nodal values by plotting them using the "caper" plot function (Felsenstein, 1985).
Subsequently, we investigated the allometric relationship between log (male) (dependent variable) and log (female) (independent variable) body size, specifically thorax length, and width, femur III length, and femur I width, by fitting four major axis regressions (model II regression; MA, Sokal & Rohlf, 1981) utilizing the phylogenetic independent contrasts (Garland et al., 1992).Rensch's rule predicts the slope of male on female size to be significantly steeper than 1.To account for the expected mean value of contrasts being zero, we enforced the MA regression to pass through the origin (Sol et al., 2014).We reported the major axis regression slopes (β), along with their corresponding 95% lower and upper confidence intervals, calculated using the "smatr" R package (Warton et al., 2006(Warton et al., , 2012)).

| Phylogenetic reconstruction
We obtained a fully resolved phylogeny of Pyrgomorphidae (Figure 3), with most phylogenetic relationships highly supported (bootstrap values ≥95%

| Ancestral state reconstruction of body size
We obtained morphological information for 32 grasshopper species that were used to perform ancestral state reconstruction of body size.
This sample includes species from Africa, America, Australia, and Asia.
Based on the current range of pyrgomorphid size variation, the body size reconstruction suggests that the ancestral pyrgomorphid likely had an intermediate size.The femur III and I were smaller than those of the present species, while the thorax length and width were larger than the average in the present species (Table 5).Overall, the group exhibited a trend toward increasing body size.However, there have been at least four independent times where a trend toward evolving a smaller body size has occurred (Figure 4a-d).Approximately, 32% of the nine species analyzed showed a trend toward the reduction of their body size: Psedna nana, Algete brunneri, Jaragua oviedensis, Piscacris robertsi, Ichthyotettix mexicanus, Sphenotettix nobilis, Pyrgotettix pueblensis, Sphenacris crassicornis, and Pyrgomorpha conica.

| Climatic conditions and body size
According to the comparisons between the pgls models with the AIC criteria, the best predictor model was only associated with the mean annual temperature, and the thorax width SSD (F (3,28) = 3.385, p = .032).The mean annual temperature showed a positive relationship with the thorax width SSD (see Table 6).Thus, lower annual temperatures tend to favor the presence of species with relatively smaller females, while environments with higher temperatures are more likely to be inhabited by species exhibiting the highest sexual size dimorphism (SSD) bias toward females.Moreover, the divergence between the sexes is not constrained by the phylogenetic history of the group (λ = .66;p = .18).

| Body size differences between winged and wingless grasshoppers
We found specimens for 10 winged and 22 wingless species.All the branches of the phylogeny are represented by winged and wingless species (Figure 3), suggesting that throughout the evolutionary history of pyrgomorphids, the loss of wings has been recurrent.There were significant body size differences between wingless and winged males and females.However, opposite to our hypothesis winged males and females were larger than males and females of wingless species (Figures 5 and 6).After controlling for phylogenetic effects, just the femur III length and the thorax length of the males remain significant (Table 7), while for the females remain significant the F I, thorax with and the femur III length (Table 8).Moreover, we found a weak, but statisctically significant phylogenetic signal for wing condition on the family (δ = .952,p-value = .02),which is consistent with the scattered distribution of the winged condition across the phylogeny with few clades mainly comprising winged species (Figure 3).

| Rensch's rule
Our results indicate strong coevolution between males and females (Table 9).Males' traits showed slopes greater than 1.0 when regressed on females' traits, with femur III and thorax lengths significantly differing from isometry (β = 1.0; Figure 7).

| DISCUSS ION
Considering the current variation in body size of the members of the family Pyrgomorphidae, our results suggest that these grasshoppers have evolved from an ancestor of intermedium size.The absence of phylogenetic constraints (λ = .66;p = .18)has not limited the diversification of the body in this linage of grasshoppers, allows for a greater flexibility in adapting to mean annual temperatures, which may contribute to the observed differences in body size between males and females.Furthermore, the disparity in wing development in pyrgomorphs is closely linked to variations in body size.This correlation suggests that the evolution of wings in these insects is intricately tied to their physical dimensions.As body size differs among pyrgomorph species, so too does the development of their wings.
This phenomenon underscores the significance of body size in shaping the morphological characteristics of these insects.The insight into the connection between body size and wing development sheds light on the adaptive strategies employed by pyrgomorphs and provides valuable information for further research into the evolutionary mechanisms governing insect morphology.
Female fecundity and male mating success tend to increase with body size, but the optimal body size often differs between males and females, generating sexual size dimorphism (SSD).The differences between females and males are further influenced by climatic conditions such as temperature and season length, which affect growth rates and development times.Seasonality plays a crucial role in determining the duration of optimal developmental conditions for both males and females.As a result, the body size differences between males and females within pyrgomorph species are not only shaped by reproductive strategies, but also by environmental factors that regulate growth and development.

| Ancestral state reconstruction of body size
The body size reconstruction suggests that the members of the family Pyrgomorphidae have evolved from an ancestor of the intermedium size.During the evolution of the groups, there have been trends to increase and decrease body size.Both trends could be modeled by several selective pressures associated with the environments that occupy the species of the family.Small body size often correlates with a high seasonal environment, where F I G U R E 3 Dated phylogeny of Pyrgomorphidae and wing condition of the included species.The phylogeny is based on maximum likelihood analysis of mitochondrial genome sequences of 35 species (32 Pyrgomorphidae and 3 outgroup species).Black circles on the tree represent nodes with bootstrap support values ≥95%, nodes without black circles had lower support (bootstrap <80%).Colored boxes at the tips indicate wing condition, red for winged and blue for wingless species.
TA B L E 5 Body size differences between the hypothetical ancestor of Pyrgomorphidae and the average size for the femur III length, femur I width, thorax length, and thorax width of the 33 species of the family measured for this study.

F I G U R E 4
Ancestral phylogenetic reconstructions of (a) femur I width, (b) femur III length, (c) thorax length, and (d) width for 32 pyrgomorphid grasshoppers' taxa performed in R package "phytools" (Revell, 2013).For the analysis, we used the ultrametric phylogeny and the values of SDI estimated for each species.
| 13 of 19 CUEVA del CASTILLO et al. tend to increase and decrease body size.Interestingly, Pyrgomorpha conica is the only small species with wings represented in the phylogeny, perhaps its wide distribution across South Europe, North Africa, and Central Asia may be attributed to its high dispersion capacity associated with flight.In any case, we must point out that the limited number of species considered in this phylogenetic reconstruction could limit the accurate interetation of the ancestral state of the lineage and its evolutionary trends.

| Climatic conditions and body size
The impact of natural and sexual selection on the body size divergence between males and females of pyrgomorphids appears to vary across different environmental conditions.Interestingly, our findings indicate that it is the mean annual temperature rather than precipitation that plays a significant role in this divergence.
We observed that species with less sexual dimorphism tend to inhabit cooler environments, whereas those with a greater bias toward females are found in warmer regions.These patterns suggest a strong relationship between maturation time and body size in these grasshoppers (Cueva del Castillo & Nuñez-Farfán, 1999, 2002).As body size is strongly correlated with maturation time in insects, faster sexual maturation is reached at the expense of having a small body size (Thornhill & Alcock, 1983).However, in Orthoptera and Coleoptera there are species in which the instar number is higher in favorable conditions.(Esperk, Tammaru, & Nylin, 2007;Esperk, Tammaru, Nylin, & Teder, 2007).Interestinlgly, Sphenarium histrio, a species represented in this study, male and female collected from low altitude sites in the field and reared in the laboratory were larger than those from a high altitude.
Grasshoppers from a high altitude hatched earlier, had a shorter development time, presented fewer instars, and were smaller than grasshoppers from the other sampling sites (Ramírez-Delgado & Cueva del Castillo, 2022).This local adaptation highlights the intricate interplay between environmental factors and the evolution of sexual dimorphism in these insects.However, in the absence of phylogenetic constraints, when fecundity selection outweighs sexual selection, it is likely that larger female body sizes will result compared with males.

| Body size differences between winged and wingless grasshoppers
Opposite to our hypotheses, winged males and females were larger than wingless species, suggesting that the costs to produce and maintain the wings do not impact the resources canalized to body size.In fact, it is possible that if there are these costs for the winged species, maybe they are compensated for an increase in the resources canalized to body size.Alternatively, given the toxicity of some of the largest species (Yang et al., 2019), lower predation levels may lead to longer life spans and larger body sizes, indirectly benefiting fecundity.This effect could be more pronounced in warmer and less seasonal environments, where natural selection favors these traits (Horne et al., 2015).

| Rensch's rule
Despite the potential tradeoffs between maturation time and body, and the size differences between winged and wingless grasshoppers, positive directional sexual selection on femur III and thorax lengths may provide insight into the evolution of SSD in Pyrgomorphidae, with male size hyperallometry reflecting their greater evolutionary divergence due to strong mate competition.
Larger males have advantages in accessing females and resisting

TA B L E 9
Mayor axis regressions of the independent contrast of male traits on the female traits, slope values (β) lower and upper confident intervals, variance explained by the models and p-values of the regression, and the deviation of β = 1 (H 0 ) are shown.& Nuñez-Farfán, 1999;Wickler & Seibt, 1985), with higher mating success and extraordinary copulation and mate-guarding durations reported in some species (Cueva del Castillo, 2003;Descamps M & Wintrebert D., 1966;Wickler & Seibt, 1985).These behaviors may decrease the possibility of sperm competition.

F
Mean and standard error values for (a) femur I width, (b) femur III length, (c) thorax width and (d) thorax length of the wingless and winged females of the family Pyrgomorphidae.F I G U R E 6 Mean and standard error values for (a) femur I width, (b) femur III length, (c) thorax width, and (d) thorax length of the wingless and winged males of the family Pyrgomorphidae.TA B L E 7 Phylogenetic analyses of variance for (a) femur III length, (b) femur I width, (c) thorax length, and (d) thorax width of wingless and winged males of the family Pyrgomorphidae.
The family Pyrgomorphidae shows great diversification, successfully colonizing a variety of niches across their distribution range from Asia to Tropical America(Sanabria & Cueva del Castillo, 2020).This broad range of habitats suggests that the relative impact of natural and sexual selection on body size has likely shifted in response to changing environmental conditions, that favored the body size divergence of females and males.Variations in pyrgomorphids may be closely linked to their evolution in distinct continents.For example, it is noteworthy that toxicity is absent in numerous neotropical species of the Genus Sphenarium and Prosphena.This suggests that environmental factors and geographical isolation have played a significant role in shaping the characteristics of these species.Understanding these differences can provide valuable insights into the intricate processes of evolution and adaptation.The specific causal mechanisms driving these changes in SSD remain an open question that requires further investigation.In any case, further investigation is required to determine the causal mechanisms behind the changes in the direction of SSD in this group of grasshoppers.AUTH O R CO NTR I B UTI O N S Raúl Cueva del Castillo: Conceptualization (equal); data curation (lead); formal analysis (equal); funding acquisition (lead); investigation (lead); methodology (lead); project administration (lead); resources (equal); software (equal); supervision (lead); validation (equal); visualization (lead); writing -original draft (lead); writingreview and editing (equal).Salomón Sanabria-Urbán: Data curation (equal); formal analysis (equal); investigation (equal); methodology (equal); software (equal); validation (equal).Ricardo Mariño-Pérez: Data curation (equal); investigation (equal); validation (equal); writing -review and editing (equal).Hojun Song: Data curation (equal); formal analysis (equal); methodology (equal); software (equal); supervision (lead); validation (lead); writing -review and editing (equal).ACK N OWLED G M ENTS R C del C is grateful to UNAM DGAPA-PASPA program for a sabbatical fellowship and CONACYT for complementary support (Program: Apoyos complementarios para Estancias Sabaticas Vinculadas a la Consolidacion de Grupos de Investigacion) at Texas A & M University.The authors are grateful for the invaluable support offered by Greg Cowper during our consultation of the Philadelphia Academy of Natural Sciences pyrgomorphid specimens, and to VH Jìmenz-Arcos for his valuable help running some of the analyses.The authors are grateful for the invaluable support offered by Greg Cowper during our consultation of the Philadelphia Academy of Natural Sciences pyrgomorphid specimens.CO N FLI C T O F I NTER E S T S TATEM ENTNone declared.DATA AVA I L A B I L I T Y S TAT E M E N TThe GenBank accession numbers of the used sequences are shown in TableS1and the final alignment used and the resulting phylogeny are available at TreeBASE under submission 31,267.

F I G U R E 7
Significant allometric major axis regressions of independent contrasts (IC) of (a) the femur III length of males on the femur III length of females', and (b) the thorax length of the males on the thorax length of the females.Dashed lines indicate isometry (β = 1).
Lovich and Gibbons (1992)sexual size dimorphism index (SSD) of the four morphological traits measured to males and females of 32 Pyrgomorphid species.
(Song et al., 2020)largelyZahid et al., 2021)evious phylogenetic studies (Mariño-Pérez & Song, 2019;Zahid et al., 2021)and phylogenomic inferences(Song et al., 2020)on this group, with minor differences in nodal support values.TA B L E 2 Taxonomic information and Genbank accession numbers for the 35 taxa used in the phylogenetic analysis.
TA B L E 6 (a) Parameters included in the significative multiple regression model of the SSDi of the thorax width of pyrgomorph grasshoppers on the climatic variables and (b) Delta Akaike Information for different regression models used to explain variation in the thorax width SSD in pyrgomorphid grasshoppers as a function of climatic variables.
Note: Results in bold represent the best regression model as indicated by ∆ AIC.SE, t and p-values are shown.Abbreviations: ∆AIC, delta Akaike information criterion; A, annual mean temperature; AIC, Akaike's information criterion; B, temperature seasonality; C, mean temperature of the wettest quarter of the year; D, annual precipitation; E, precipitation seasonality; F, precipitation of the wettest quarter of the year; r 2 , variance explained by each model; SE, standard error.
Phylogenetic analyses of variance for (a) femur III length, (b) femur I width, (c) thorax length, and (d) thorax width of wingless and winged females of the family Pyrgomorphidae.
Abbreviations: df, degrees of freedom; F, F values; MS, mean squares, p, p-value; SS, sum of squares.*p-valuecorrected by the phylogeny (p*) are shown.TA B L E 8Abbreviations: df, degrees of freedom; F, F values; MS, mean squares; p, p-value; SS, sum of squares.*p-valuecorrected by the phylogeny (p*) are shown.